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Abstract. A 3-dimensional model of the synchrotron 
emissivity of the cocoons of powerful extragalactic radio 
sources of type FRII is constructed. Bulk backflow and 
energy losses of the relativistic electrons, radiative and 
adiabatic, are self-consistently taken into account. Thus 
the model is an extension of spectral aging methods in- 
cluding the underlying source dynamics into the age es- 
timates. The discrepancies between spectral ages and dy- 
namical ages arising from earlier methods are resolved. It 
is also shown that diffusion of relativistic particles within 
the cocoon is unlikely to significantly change the particle 
spectrum and thus the emitted radio spectrum. Projec- 
tion along the line of sight yields 1 or 2-dimensional sur- 
face brightness distributions which can be compared with 
observations. From the model parameters constraints on 
the source age, the density of the source environment and 
the angle to the line of sight can be derived. Application 
of the method to Cygnus A, 3C 219 and 3C 215 show that 
the method provides robust estimates for the model pa- 
rameters for sources with comparatively regularly shaped 
radio lobes. The resolution of the radio maps required is 
only moderate. Within the large uncertainties for the ori- 
entation angle, the three example sources are found to be 
consistent with orientation-based unification schemes of 
radio- loud AGN. In the case of Cygnus A the gas density 
of the environment is found to agree with independent 
X-ray measurements. For 3C 219 and 3C 215 the densi- 
ties derived from the model are apparently too low. It is 
suggested that these discrepancies are caused by overes- 
timates of slope and core radius of /3-models for the gas 
density distribution from X-ray observations for clusters 
hosting powerful radio sources. 

Key words: Galaxies: active - Galaxies: individual: 
Cygnus A - Galaxies: individual: 3C 219 - Galaxies: jets - 
Quasars: individual: 3C 215 - Radio continuum: galaxies 
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1. Introduction 

The 'standard model' of extragalactic classical double 
radio sources by Scheuer (1974) explains these sources 
as twin jets emerging from the Active Galactic Nucleus 
(AGN) and impinging on the surrounding gas. The com- 
pact, high surface brightness regions or hot spots at the 
end of the jets are interpreted as the sites of the interac- 
tion between the jets and the environment. After passing 
through the hot spot region, the jet material inflates the 
radio lobes or cocoon which is observed as diffuse emission 
in between the hot spots and the source core (e.g. Muxlow 
& Garrington 1991 and references therein). This picture 
forms the basis of virtually all more recent attempts at 
modeling the dynamics and radio emission properties of 
powerful radio galaxies and radio-loud quasars (Begelman 
& Cioffi 1989, Falle 1991, Daly 1994, Nath 1995, Kaiser & 
Alexander 1997, Chyzy 1997, Kaiser et al. 1997, Blundell 
et al. 1999). 

These models for the evolution of individual classical 
doubles or FRII-type objects (Fanaroff & Riley 1974) can 
be used to study the cosmological evolution of the pop- 
ulation as a whole. One of the more important trends is 
the apparent decrease of the mean linear size of the radio 
structures with increasing redshift (see Neeser et al. 1995 
and references therein). 

Two recent attempts in fully explaining the linear size 
- redshift anti-correlation were presented by Blundell et 
al. (1999) and Kaiser & Alexander (1999a). Blundell et al. 
argue that a specific form of pre-aging of the relativistic 
particle population in the hot spots in connection with the 
lower flux limit of complete observed samples causes this 
anti-correlation (the so-called youth-redshift degeneracy, 
Blundell et al. 1999, Blundell & Rawlings 1999). Alter- 
natively, Kaiser & Alexander propose that the apparent 
smaller sizes of radio sources at high redshift could be 
caused by a denser environment of these objects at high 
redshift. Clearly, the evolution or non-evolution of the ra- 
dio source environments is of great interest to decide this 
and other important cosmological questions. 

Various methods have been employed in determining 
the properties of the radio source environments. X-ray ob- 
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servations of the hot gas around radio sources in clusters 
yield direct estimates of the gas density (Crawford et al. 
1999, Hardcastle & Worrall 1999). However, such studies 
are at present confined to objects at low redshifts. Fur- 
thermore, the AGN (e.g. Sambruna et al. 1999) and the 
large-scale radio structure contribute to the X-ray emis- 
sion (Brunetti et al. 1999, Kaiser & Alexander 1999b, Sect. 
^) . The properties of the density distribution inferred from 
X-ray observations may therefore be considerably influ- 
enced by the presence of the radio source itself. 

Faraday rotation of the polarisation angle of the syn- 
chrotron emission and the related depolarisation can be 
used to determine the gas density of the material sur- 
rounding the radio lobes (e.g. Garrington et al. 1988, 
Laing 1988). Unfortunately, this method does not provide 
a direct measure of the gas density as the rotation mea- 
sure also depends on the strength of the magnetic field in 
the source environment. Usually it is not possible to break 
this degeneracy because the strength of the magnetic field 
is only poorly known. 

Constraints on the density of the radio source environ- 
ments also come from optical or infrared galaxy counts 
around the host galaxies (e.g. Hill & Lilly 1991). It is not 
straightforward to decide whether galaxies in the field of 
the radio source host are associated or chance background 
objects. Resolving the ambiguity would ideally require the 
spectroscopic measurement of the redshifts of all objects 
in question. This is very time consuming. In any case, the 
method provides only indirect constraints on the gas den- 
sity around the radio source as this has to be inferred from 
comparison with low redshift clusters or groups of similar 
richness. 

The ages of powerful radio sources can in principle be 
determined from their radio spectrum (e.g. Alexander & 
Leahy 1987). The various energy losses of the relativistic 
particles depend on time and so the shape of the radio 
spectrum contains an encoded history of the source. In 
practice the time-dependence of the energy losses compli- 
cates the estimation of the spectral age because different 
parts of the source have different ages. Even when radio 
maps at various frequencies which fully resolve the radio 
lobes are used, it is difficult to disentangle the effects of 
the various loss mechanisms and possible bulk backflow 
of the cocoon material along the jet (e.g. Rudnick et al. 
1994). 

The model developed in this paper aims at tracing the 
individual evolution of parts of the cocoon and thereby 
providing more accurate estimates for the source age. At 
the same time the model also constrains the density in the 
source environment and other parameters like the energy 
transport rate of the jets. It is solely based on radio obser- 
vations which are available for a large number of objects, 
even at high redshift. This model may therefore provide an 
important step in determining the cosmological evolution 
of the FRII radio source population. 



In Sect. [2] I show that diffusion of relativistic particles 
in the cocoons of FRII-type objects should not signifi- 
cantly change the energy distribution of particles. In Sect. 
U the dynamical model of Kaiser & Alexander (1997, here- 
after KA) and its extension by Kaiser et al. (1997, here- 
after KDA) to include synchrotron emission are briefly 
summarised. A 3-dimensional model of the synchrotron 
emissivity of the cocoon based on this analysis is con- 



structed in Sects. 3.3 and 3.4, Methods for comparing the 



model predictions with observations are developed in Sect. 
[|. The degeneracy of model parameters resulting from the 
comparison method is also discussed here. The model is 
then applied to three FRII-type radio sources, Cygnus A, 
3C 219 and 3C 215, in Sect. |. The results are discussed 
in Sect. [| The main conclusions are summarised in Sect. 

I- 

Throughout this paper I use H a — 50 km s _1 Mpc -1 
and q Q = 0. 



2. Comparison to spectral aging methods 

The age of extragalactic radio sources can be estimated 
by the use of spectral aging arguments. This method re- 
lies on the determination of the break in the radio spec- 
trum caused by the time-dependent energy losses of the 
relativistic electrons within the cocoon (e.g. Alexander & 
Leahy 1987). At frequencies higher than the break fre- 
quency, 1^,, the spectrum significantly steepens due to the 
radiative energy losses of the electrons. The spectral age of 
a electron population, i sa , with observed break frequency 
i>b within a magnetic field of strength B which does not 
vary in time is proportional to B~ z l 2 v h 1 . Despite the 
implication of source dynamics that the magnetic field 
does vary in time, the constant field relation is usually 
used in determining the spectral age of a given source. 
In general the strength of the magnetic field will decrease 
while the source expands and so t sa will be an overestimate 
of the true source age, t. However, because of variations of 
the magnetic field in the cocoon and the backflow of gas 
in this region, t sa is usually found to be lower than the 
dynamical ages inferred from the advance speed of the 
cocoon by various methods (Alexander & Leahy 1987). 

The model presented in the following can be viewed as 
an extension of the spectral aging formalism. Within the 
cocoon it traces the evolution of the local magnetic field in 
time. This allows the accurate determination of the energy 
distribution function of the relativistic particle population 
and thus the emitted spectrum at a given location along 
the cocoon. It is therefore not surprising that the model 
predicts an older age for Cygnus A (see Sect. |j) than the 
classical spectral index analysis of Carilli et al. (1991). 
The combination of a dynamical model with the accurate 
treatment of the local evolution within the cocoon of the 
radio emission properties implies that the spectral and 
dynamical ages are identical. 
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.1. Diffusion of relativistic particles 



3.1. The dynamical model 



In the context of spectral aging methods other processes 
changing the energy distribution of the relativistic elec- 
trons and thus invalidating the age estimates have been 
put forward. These are summarised and discussed by 
Blundell & Rawlings (2000). They show that most of these 
processes are rather inefficient and will not strongly in- 
fluence spectral aging methods or the model discussed 
here. However, Blundell & Rawlings (2000) claim that 
the anomalous diffusion mechanism of Rechester & Rosen- 
bluth (1978) can lead to very fast diffusion of relativistic 
particles through the tangled magnetic field in the cocoon. 
So much so that the radio spectra observed at any point 
along the cocoon essentially arise from the electrons accel- 
erated by the jet shock at the hot spot within the last 10 6 
years of the observation. This would imply that no infor- 
mation on the source age can be derived from the spatial 
properties of the observed emission. Support for a univer- 
sal particle energy spectrum may come from the interpre- 
tation of observations of Cygnus A by Katz-Stone et al. 

(1993) , Rudnick et al. (1994) and Katz-Stone & Rudnick 

(1994) . They find that a single, non-power law spectrum 
shifted in frequency by the local strength of the magnetic 
field is emitted by all parts of the lobes of Cygnus A. Note 
however, that the spectral shape they find may be caused 
by the free-free absorption in our own galaxy of the emis- 
sion of Cygnus A at low frequencies (Carilli et al. 1991. In 
this case, the universal spectrum is simply explained by 
shifting an aging energy spectrum of relativistic particles 
which is not a single power law. Alternatively, changes in 
the strength of the magnetic field along a given line of sight 
may also introduce additional curvature in the observed 
spectrum (Rudnick et al. 1994). The mentioned interpre- 
tation of the radio observations of Cygnus A in itself is 
therefore no proof of efficient diffusion acting in the lobes 
of radio galaxies. 

In the appendix I show that anomalous diffusion is 
probably much less effective in the cocoon plasma as pre- 
viously thought. In general we do not observe any signa- 
ture for diffusion losses of the cocoons of radio sources 
and so diffusion will not alter the distribution of relativis- 
tic particles within the cocoon. This allows us to use the 
spatial distribution of the synchrotron radio emission of 
FRII sources to infer their age. The model developed in 
the following can be viewed as an extension to the clas- 
sical spectral aging methods in that it takes into account 
the evolution of the magnetic field in the lobe. 



3. The model 

In this Sect. I briefly summarise the dynamical and radia- 
tive model that form the basis for the extended treatment 
presented in this paper. Following this the prescription for 
the spatial distribution of the synchrotron emission within 
the radio lobes is developed. 



The large scale structure of radio galaxies and radio-loud 
quasars of type FRII is formed by twin jets emanating 
from the central AGN buried inside the nucleus of the 
host galaxy. The jets propagate into opposite directions 
from the core of the source. They end in strong jet shocks 
and, after passing through these shock, the jet material 
inflates the cocoon surrounding the jets. The cocoon is 
overpressured with respect to the ambient medium and 
therefore drives strong bow shocks into this material. 

Falle (1991) and KA showed that the expansion of the 
bow shock and the cocoon should be self-similar which is 
supported by observations (e.g. Leahy & Williams 1984, 
Leahy et al. 1989, Black 1992). In these models the density 
distribution of the material the radio source is expanding 
into is approximated by a power law, p x = p (r /a )~@ , 
where r is the radial distance from the source centre and 
a Q is the core radius of the density distribution. X-ray 
observations of groups and clusters of galaxies show that 
the density of the hot gas in these structures is often dis- 
tributed according to (e.g. Sarazin 1988) 
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1 + (r/a f 



3/3'/2 ' 
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Outside a few core radii, a , the power law assumed above 
with (3 — 3/3' provides an adequate fit to Eq. (|l]) . Even for 
smaller distances r good power law approximations can 
be found by adjusting j3 and p (e.g. Kaiser & Alexander 
1999a). 

In the model of KA it is also assumed that the rate at 
which energy is transported along each jet, Q , is constant 
and that the jets are in pressure equilibrium with their 
own cocoon. The very high sound speed within the cocoon 
results in a practically uniform pressure within this region 
apart from the tip of the cocoon. The pressure in this 
'hot spot' region, named for the very strong radio emission 
originating in the shock at the end of the jets, is somewhat 
higher as the cocoon material injected by the jets at these 
points is not yet in pressure equilibrium with the rest of 
the cocoon. 

In the following I will concentrate on only one jet and 
the half of the cocoon it is contained in. From KA I take 
the expressions for the evolution of the uniform cocoon 
pressure, 
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Here, T x is the ratio of specific heats of the gas surrounding 
the radio source, t is the age of the jet flow and C\ is a 
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dimensionless constant. The ratio, Ph c , of the pressure in 
the hot spot region and p c is constant in the model of KA. 

KDA extend the model of KA to include the syn- 
chrotron emission of the cocoon. This is done by splitting 
up the cocoon into small volume elements, SV C , the evo- 
lution of which is then followed individually. By assuming 
that these elements are injected by the jet into the hot 
spot region at time ti during a short time interval Sti and 
then become part of the cocoon, KDA find 



SVc 



(To - 1) Qo 

Pc (ti) 



(P hc ) 



(i-r c )/r c ft 



(4) 



where T c is the ratio of specific heats for the cocoon mate- 
rial and C4 = (4 + (3) / [T c (5 — f3)]. Because the expansion 
of the cocoon is self-similar we can set for the volume of 
the total volume of the cocoon, following the notation of 
KA, V c = C3L?, where C3 is a dimensionless constant and 
depends on the geometric shape of the cocoon. In order to 
ensure self-consistency the integration of Eq. (0) over the 
injection time from ti = to ti = t must be equal to the 
total cocoon volume V c . Substituting Eq. (||) for Lj then 
yields 
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Note, that this expression for c\ is different from the one 
given by KA. In their analysis they used the expression 
for conservation of energy for the entire cocoon 



dU c 
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[d( Pc V c ) + d( Ph V h )} = Q Q dt, 



(6) 



where 14 is the volume of the hot spot region with pres- 
sure p^. KA then used the simplifying assumption that the 
cocoon has a cylindrical geometry, the expansion of which 
is governed along the jet axis by ph while its growth per- 
pendicular to this direction is driven by p c . This implies 
Ph/ P c = Phc = 4i?|., where Rt is the ratio of the length of 
one jet and the full width of the associated lobe halfway 
down the jet. Kaiser & Alexander (1999b) subsequently 
derived empirical fitting formulae for Ph c as functions of 
/3 and Pt from an analysis of the flow of shocked gas be- 
tween the bow shock and the cocoon. Their results showed 
that the original approximation tends to overestimate the 
value of Phc • In the following I use a generalised empirical 
fitting formula based on their result and further calcula- 
tions with additional values of 0, 



P hc = (2.14-0.52/3) R T - 



2.04-0.25/3 



(7) 



In order to satisfy Eqs. (||) and (jfj) I now generalise the 
approach of KA by setting d(p\ l V\ l ) — fd(p c V c ). Because 
of the self-similar expansion of the cocoon / is a constant 
and will depend on the ratio Ph c - The strict separation of 
the cocoon material into the hot spot region and the rest 
of the cocoon with two distinctly different values of the 



respective pressure within these regions is of course artifi- 
cial. In the cocoons of real FRII objects the transition of 
jet material from hot spot to cocoon will be accomplished 
in a continuous hydrodynamical flow along a pressure gra- 
dient much smoother than the sudden change from ph to 
p c described here. However, a detailed model of this flow 
is beyond the scope of this paper and for simplicity I will 
use the assumption of a strict spatial separation in the fol- 
lowing. As the two cocoon regions are in physical contact 
with each other, / may become negative as the expansion 
of one region may influence the evolution of the other. 
Substituting Eqs. (|j) and (||) into Eq. (|h gives 
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3.2. Synchrotron emission 

The cocoon volume elements SV C are filled with a mag- 
netised plasma and a population of relativistic electrons 
accelerated at the shock terminating the jet flow at the 
hot spot. They therefore emit synchrotron radio radiation. 
In optically thin conditions the monochromatic luminos- 
ity due to this process can be calculated by folding the 
emissivity of single electrons with their energy distribu- 
tion (e.g. Shu 1991). 

Following KDA I assume that the initial energy dis- 
tribution of the relativistic electrons as they leave the ac- 
celeration region of the hot spot follows a power law with 
exponent — p between 7 = 1 and 7 = 7 max . The electrons 
are subject to energy losses due to the adiabatic expansion 
of SV C , the emission of synchrotron radiation and inverse 
Compton scattering of the CMB. For a given volume ele- 
ment at time t that was injected into the cocoon at time 
ti these losses result in an energy distribution (see KDA) 

-p/2 / , \ c 4 (p+2)/3 
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where u cm b (z) is the energy density of the CMB radiation 
field at the source redshift z and m c is the rest mass of an 
electron. Here I have used cq = 1 — C4 (T b + 1/3) and C7 = 
1 — C4/3. The normalisation of the energy spectrum, n (t{), 
is given by integrating the initial power law distribution 
over the entire energy range 

n (ti) = j 
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or, for p = 2, 

U e (ti) 



n (ti) 



log (7max) + 



(11) 



where u c is the total energy density of the relativistic par- 
ticle distribution. At time ij I set ue(ti) / 'u e {t\) — r and 
for simplicity p c {h) = (r c — 1)[ub(^i) + u e {t{)]. In the 
following I assume that the minimum energy condition 
(e.g. Miley 1980) is initially fulfilled in each volume ele- 
ment SV C and therefore r = (p + l)/4. From Eqs. (§) and 



it follows that p c (tj) = p c (t) (t/^i 



c 4 r c 



and with the 



assumption of completely tangled magnetic fields I find 
ub (t) = u B (U) (</ti) C4rB (see also KA). With this the set 
of equations describing the radio synchrotron emissivity, 
— Pv/6Vc, of a given volume element injected into the 
cocoon at time t[ only depends on the present value of 
the pressure in the cocoon, p c (t), and the age of the radio 
source, t. 

3.3. Spatial distribution of the emission 

So far the cocoon volume elements SV C were only charac- 
terised by their injection time into the cocoon, t\. Prom 
the analysis above it is not possible to decide where they 
are located spatially in the cocoon. 

From the above analysis it is clear that the radio spec- 
trum emitted by a given cocoon volume element depends 
on the 'energy loss history' of this part of the cocoon. Adi- 
abatic losses only change the normalisation of the emitted 
spectrum while its slope at a given frequency is governed 
by the radiative loss processes of the relativistic electrons. 
In the model described above the strength of the mag- 
netic field which determines the magnitude of synchrotron 
losses is tied to the pressure in the cocoon, p c . The value 
of p c in turn depends on the energy transport rate of the 
jet, Q , and a combination of parameters describing the 
density distribution of the gas the cocoon is expanding 
into, p cio- In the analytical scenario presented here the 
volume elements are the building blocks of the cocoon and 
the variation of the radio surface brightness of the cocoon 
of an FRII along its major axis can therefore potentially 
provide information on the properties of the source envi- 
ronment. For this I now identify the cocoon volume ele- 
ments SV C with infinitesimally thin cylindrical slices with 
their radius, r, perpendicular to the jet axis. A similar 
approach is used in the model of Chyzy (1997). However, 
in this model the dynamical evolution and the radio lu- 
minosity of the entire cocoon depend on the geometrical 
shape and evolution of the volume elements. In the model 
of KA and KDA this is not the case and the identification 
of the 6V C with a specific geometrical shape does not alter 
their results. 

The volume of a thin cylindrical slice of the cocoon is 
given by 5V C = nr^Sz, where Sz is the very small thickness 
of the slice along the jet axis. The radius of the slices, r c , 



Table 1. Typical value ranges for the dimensionless con- 
stants in the dynamical model. This assumes < (3 < 2 
and 1.3 < R T < 6 (e.g. Leahy & Williams 1984). 



Cl 



f'3 



C4 



t'li 



C7 



[1.6,6.0] [0.01,0.3] [0.6,1.5] [-1.5,0.0] [0.5,0.; 



depends on their position along the jet axis, I, which I 
define in units of Lu the length of one half of the entire 



cocoon (see Sect. 3T). The outer edges of the slices form 
the cocoon boundary or contact discontinuity. Most FRII 
sources have cocoons of a relatively undistorted, ellipsoidal 
shape (e.g. Leahy et al. 1989). I therefore parameterise the 
cocoon boundary as 



a (l-/ Qi r, 



(12) 



where a Q , ct\ and ct2 can be determined from radio obser- 
vations of the cocoon. 

Following the observational results of Leahy & 
Williams (1984) and Leahy et al. (1989) KA and KDA 
used the aspect ratio, Rt, to characterise the geometrical 
shape of the cocoons of FRII sources. This ratio is defined 
as the length of one side of the cocoon measured from the 
radio core to the cocoon tip divided by its width mea- 
sured half-way along this line. Using this definition it is 
straightforward to express a Q in terms of Rt as 



Qo 



2i? T 



(13) 



The dimensionless volume constant C3 now becomes 



a 



2 rl 



C3 



w 7? / (l-^ 1 ) 2 " 2 dl 
L i Jo 

a 2 

7T^2 B ( 2 «2 + l;l/«l) 



(14) 



where B(/i; v) is the complete Beta-function. 

Tab. |l| shows typical values for the dimensionless con- 
stants in the model. 



3.4. Backflow 

The cylindrical slices are injected into the cocoon at a time 
U. For simplicity I assume that the slices remain and thus 
move within the cocoon as entities afterwards. In other 
words, I neglect any mixing of material between slices and 
I also assume that the geometrical shape of the slices does 
not deviate from the initial thin cylinders. Numerical sim- 
ulations (e.g. Falle 1994) show that the gas flow in the 
cocoon is rather turbulent, at least close to the hot spot 
region. It is likely that large-scale turbulent mixing in the 
cocoon leads to large distortions of the projected cocoon 
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shape as seen in radio observations. In this case, the reg- 
ular cocoon shape described by Eq. ( |l2| ) will be a poor 
representation of the 'true' cocoon shape. For such dis- 
torted sources it is unlikely that the model presented here 
will provide a good description. However, the simple pic- 
ture of cylindrical slices may still represent the 'average' 
behaviour of the gas flow in more regularly shaped cocoons 
rather well. 

The model is designed to constrain source and envi- 
ronment parameters using mainly the gradient of the ra- 
dio surface brightness along the cocoon. Problems with 
this simplified model will therefore arise if the relativistic 
electrons in the cocoon are distributed efficiently by dif- 
fusion. In Sect. H I show that the diffusion of relativistic 
particles is unlikely to change their distribution on large 
scales. It is therefore reasonable to assume that the rel- 
ativistic particles are effectively tied to the cocoon slice 
they were originally injected into. 

Numerical simulations of the large scale structure of 
FRII sources strongly suggest that a backflow of material 
along the jet axis is established within the cocoon (e.g. 
Norman et al. 1982). The model describing the source dy- 
namics predicts the growth of the cocoon to be self-similar 
and therefore the backflow within the cocoon should be 
self-similar as well. This suggests that the position of a 
slice of cocoon material injected into the cocoon at time 
tj is given by I — (ti/t)" 3 , where a 3 governs the velocity 
of the backflow at a given position along the cocoon. In 
order for the model to be self-consistent, all SV C have to 
add up to the total volume of the cocoon, V c . Using Eqs. 
(|), (|), d), © and @ and replacing U/t by l 1 ^ 3 = x 
an implicit Eq. for a?, can be found from this integration; 



air c (5-/3) 
(9r c -4-/3) 



B (2a 2 + 1; i/«i) 1 
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> (1 - x aia3 )~ 2a2 dx 
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aia.3 



-B (1 - 2a 2 ; [c 4 (T c - 1) + 1] / (c*ia 3 )) 



Note that this expression requires < a 2 < 1/4 and gives 
«3 > 0. For the values of the shape parameters used in 
Sect. | (on = 2 and a 2 = 1/3) and /3 = 1.5 I find a 3 ~ 4.7. 
Note that 0-3 — > 00 for a 2 — > 1/4. 

The backflow velocity within the cocoon is given by 



fback = t I 1 — a 3 I Lj, 



(16) 



where a dot denotes a time derivative. In the rest frame 
of the host galaxy the backflow is observed to flow in the 
direction of the source core for tWk < 0. For Wback = the 
cocoon material is stationary in this frame and for positive 
values the backflow is strictly speaking not a 'backflow' 
but trailing after the advancing hot spot. 



The backflow is fastest just behind the hot spot and 
decelerates along the cocoon. The deceleration implies a 
pressure gradient along the cocoon which may seriously 
violate the assumption made for the dynamical model of 
a constant pressure throughout the cocoon away from the 
hot spots. To estimate the magnitude of the pressure gra- 
dient I use Euler's equation 



dv 



back 



1 dp (I) 



dt 



Pc 



Ljdl 



(17) 



The density of the cocoon material, p c , is given by 



QoSti 



(7j 



l)cW c 

Pc (*i) 



( 7j -i)(r c -i) c 



p(r c — i)/r c ja/013 

2 Mic 1 



(18) 



where I assumed that the entire energy in the jet is trans- 
ported in the form of kinetic energy of the flow with a bulk 
velocity corresponding to the Lorentz factor 7j. With this, 
it is straightforward with the help of Eq. ( |l6| ) to integrate 
Eq. ( |l7| ) which yields the pressure along the cocoon as a 
function of I: 



P(l) 



= 1 - 



(a' 3 -l)(3a' 3 + 2-/3)P 1 



he 



-i)/r„ 



L 



3( 7j -i)(r c -i)[2 + C4 A*3(i-r c )] 
\ 2 



J ] ^2 + c 4 /Q3(i-r c ) 



(19) 



where a' 3 



(5 — (i) a 3 /3 and I have used the condition 
p(0) — p c . The mean advance speed of the cocoon of 
an FRII source, Lj, is inferred from observations of lobe 
asymmetries to be in the range from 0.05 c (Scheuer 1995) 
to 0.1 c (Arshakian & Longair 2000). This is in good agree- 
ment with the predictions of the dynamical model used 
here (see KA). The ratio Phc is usually of order 5 (see Eq. 
0) and so even for only mildly relativistic bulk flow in the 
jet (7j ~ 2) and large gradients in the backflow velocity, 
e.g. a 3 - 10, Eq. @ predicts 0.28 <p(l) /p c < 1.0. Note 
that for jj — 5 the lower limit of this ratio increases to 
0.82. 

From this I conclude that the existence of a pressure 
gradient along the jet within the cocoon is required to 
decelerate the backflow of the cocoon material. However, 
the pressure varies only by a factor of a few at most along 
the entire length of the cocoon. This implies that within 
this limit the model is self-consistent. 

4. Comparison with observations 

The model presented here depends on a large number 
of parameters. I assume that each volume element only 
contains the relativistic particle population and therefore 
T c = T B = 4/3. The other model parameters can be 
roughly grouped into 
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— geometrical parameters: the constants describing the 
cocoon shape ai, a-i and Rt and the orientation of 
the jet axis with respect to the Line-Of-Sight (LOS), 
0, 

— properties of the initial energy distribution of the rel- 
ativistic electrons and/or positrons: the slope of the 
distribution, p, and the high energy cut-off, j max , 

— properties of the source and its environment: the pres- 
sure within the cocoon, p c , the age of the source, t, 
and the slope of the power law density distribution of 
the environment, (3. 

In the following I discuss the methods employed in com- 
paring the model predictions with observations. Using 
these methods, I will then show that because of the nature 
of the model several degeneracies between model parame- 
ters exist. Eliminating these requires further assumptions 
to be made but also reduces the complexity of the param- 
eter estimation. 



4-1- 2- dimensional comparison 

A radio map of an FRII source is composed of pixels which 
contain information on the monochromatic radio surface 
brightness of the source, S v yj) , at a given position 
yj) projected onto the plane of the sky. Here I use the 
projected distance along the jet axis, V , and perpendicular 
to the jet, y'. Both are measured in units of the projected 
length of the cocoon, In optically thin conditions gen- 
erally present in radio source lobes this projection corre- 
sponds to a LOS integral of the synchrotron emissivity, e„, 
through the 3-dimensional sour ce a t eac h p ixel location. 



Following the analysis in Sects. |3jJ and |3.2| , e v is only a 
function of the unprojected dimensionless distance from 
the core measured along the jet axis, I, for a given set of 
the eight source and environment parameters. Using the 
model described above, it is therefore possible to construct 
'virtual radio maps' for a given set of model parameters 
by projecting the 3-dimensional model onto the plane of 
the sky. Since a radio source is in general viewed with its 
jet axis at an angle 9 to the LOS, the value of I, and there- 
fore that of e u , changes along the path of the LOS integral. 
Also, the distance from the core to the tip of the cocoon 
measured on an observed radio map, Lj, is not equal to 
the physical size of the cocoon, Lj = L'-/ sin#. This, of 
course, also implies that the aspect ratio of the cocoon, 
R' T , measured in the observed map is smaller than the 
'real' value Rt = R' T / siri9, i.e. radio sources appear 'fat- 
ter' than they really are. All these projection effects have 
to be taken into account when comparing virtual radio 
maps resulting from these models to observations. 

The surface brightness is calculated at a frequency v 
which is given by the observing frequency in the frame of 
the observer, v' \ and the redshift of the source, i.e. v = (1 + 
z)v' . Finally, to account for cosmological effects the result 



of the LOS integration must be multiplied by (1 + z)/D£, 
where Dl is the luminosity distance to the source. 

The model map can be calculated at arbitrary resolu- 
tion. However, before comparing the result with the ob- 
served map it must be convolved with the beam of the 
radio telescope used for the observations. This was done 
assuming a 2-dimensional Gaussian shape for the telescope 
beam. 

Once a virtual map is compiled for a set of model pa- 
rameters this map can be compared pixel by pixel to a 
map resulting from observations of an FRII radio source 
using a \ 2 -technique, 



EE 



(20) 



Here, S v is the measured monochromatic surface bright- 
ness with rms error a v and M u is the model prediction. 

It is then possible to find the best-fitting model by 
varying the model parameters and thereby minimising the 
resulting x 2 -difference between virtual and observed map. 
The minimisation routine uses a n-dimensional downhill 
simplex method (Press et al. 1992), where n is the number 
of model parameters to be fitted. The minimisation can be 
done separately for the two halves of each cocoon since the 
model describes one jet and the associated half of the co- 
coon. Although in principle the minimisation can be done 
with one observed map at a single observing frequency, 
the constraints on the model parameters are improved by 
using two maps at two different frequencies. In this case 
the x 2 -differences resulting from the two maps are simply 
added together. In principle this is equivalent to compiling 
spectral index maps from two observed maps and compar- 
ing these with the model predictions. However, using the 
two maps directly has the advantage that pixels below the 
rms limit in one map but not in the other are not entirely 
lost for the fitting procedure. Furthermore, information on 
the absolute surface brightness in a given location is not 
contained in a spectral index map. The model would have 
to be normalised 'by hand' and the various possibilities 
to do this would lead to ambiguities in the estimation of 
model parameters. In Sec. || I use two individual maps at 
two frequency for each source. 

Using the best-fitting model parameters p c and t, the 
energy transport rate of the jet, Q , and the parameter 
combination p »o describing the density distribution in 
the source environment can be calculated from Eqs. ((H) 
and (|). 

Preliminary results using this comparison technique 
applied to radio observations of Cygnus A was presented 
in Kaiser (2000). In Sect. || this analysis is extended to 
also include 3C 219 and 3C 215. 
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4-2. 1- dimensional comparison 

The 2-dimensional comparison method described in the 
previous section requires a ray-tracing algorithm for the 
projection of the 3-dimensional model. For many maps of 
radio sources the number of pixels are so large that this 
method can become computationally very expensive. Fur- 
thermore, for maps of lower resolution the cocoon may not 
be resolved in the direction perpendicular to the jet axis. 
In such maps the surface brightness gradient along the 
jet axis can be extracted by taking a cut through the map 
along the line connecting the source core and the radio hot 
spot in the cocoon on one side. This yields a 1-dimensional 
curve of surface brightness as a function of projected dis- 
tance from the source core, S v (V). The off-axis pixels will 
not add much information. The function S v (/') can then 
be compared with the model predictions and a best-fitting 
model may be found using two maps at two observing 
frequencies as outlined in the 2-dimensional case. This 1- 
dimensional comparison involves a much smaller number 
of pixels for which a model prediction must be calculated 
than the 2-dimensional method. 



4-3. Model parameters 

To study the influence of the individual model parameters 
on the surface brightness distribution, I define a fiducial 
model with a set of fixed model parameters given in Tab. 
@. These model parameters imply a magnetic field just 
behind the hot spot of 13 nT or 130 /^G. The highest fre- 
quency of the emitted synchrotron spectrum is then ~ 300 
GHz which is close to the usually assumed high frequency 
cut-off in the spectrum of radio galaxies in minimum en- 
ergy arguments (e.g. Miley 1980). The density parameter 
p D af = 8.8 • 10 7 kg m -1 ' 5 which corresponds to a central 




0.4 0.6 
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Fig. 1. The influence of the cocoon pressure on the radio 
surface brightness along the jet axis. All curves are plotted 
for the parameters of the fiducial model except for a vari- 
ation in the cocoon pressure. The relative distance from 
the core of the source is given in units of Lj = 86.6 kpc 
(see text). The surface brightness is plotted as it would 
be extracted by an observer from a map of a source at 
redshift z = 0.1. 

density of 5 • 10~ 22 kg m~ 3 or 0.3 particles per cm -3 if 
a Q = 10 kpc. I also assume that the jet of the sources is 
100 kpc long and is viewed at an angle of 9 = 60° to the 
line of sight. This implies that it would be observed to 
have a length of Lj = 86.6 kpc corresponding to 31.2" at 
a redshift z = 0.1. For a measured R' T = 2 the aspect ratio 
of the cocoon is Ri — 2.3 for the assumed viewing angle. 

For simplicity and ease of comparison I use the 1- 
dimensional model which only predicts the surface bright- 
ness distribution along the jet axis (see previous Sect.). 
The variation of the model predictions with varying pa- 
rameters for the 2-dimensional case are essentially similar 
but the differences between maps are more difficult to vi- 
sualise. 

The solid line in Fig. |l| shows S v (I) at v' = 1.5 GHz 
assuming a pixel size of 0.3" xO. 3", corresponding to 0.8 
kpcx0.8 kpc for z = 0.1, appropriate for a telescope beam 
of 1.2" FWHM. For simplicity and in contrast to the 2- 
dimensional comparison method the averaging effects of 
the observing beam which extends over four pixels on the 
curve was not taken into account. In any case, in this Sect, 
the only interest is in gross trends of the model predictions 
when the model parameters are varied and the effective 
smoothing of the beam on the already rather smooth curve 
is small. A continuous curve is shown, since the pixel size 
is small compared to the scale of the plot. The length 
of the cocoon corresponds to more than 100 pixels. As 



pointed out in Sect. 3.2, the emission of the hot spot is 
not modeled in the approach presented here. The cocoon 
surface brightness shown in this and the following figures 
is caused by the cocoon material only after it has passed 
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Fig. 2. The influence of the source age on the radio surface 
brightness. All model parameters as in Fig. [l] except for a 
variation in the source age. 



through the hot spot region. The emission from the hot 
spot will in general dominate the total emission from the 
end of the cocoon at I — > 1. In a comparison with observed 
maps the predictions of the model can therefore not be 
used in this region. 

Also shown in Fig. [j] are the results for the same model 
with a higher and lower cocoon pressure. In the model the 
strength of the magnetic field and the energy density of 
the relativistic particles in the cocoon is tied to the co- 
coon pressure. A higher pressure therefore causes a higher 
peak of the surface brightness distribution. However, the 
increased synchrotron energy losses of the relativistic elec- 
trons also lead to a stronger gradient of the distribution 
towards the core of the source, i.e. the older parts of the 
cocoon. 

In Fig. § the effects of the source age on the distri- 
bution of S v (I) is shown. The population of relativistic 
particles at a given relative distance, I, from the core of 
the source has spent more time in the cocoon in an old 
source compared to a young source. This implies stronger 
energy losses and the peak of the surface brightness distri- 
bution is therefore located closer to the hot spot in an old 
source. Furthermore, in an observation with a given de- 
tectable threshold of S v (I) a larger fraction of the cocoon 
of a young source will be visible than of an older source. 

The same source viewed at different angles of the jet 
axis to the LOS, as shown in Fig. ||, results in a scaling 
of the surface brightness similar to the effects of chang- 
ing the cocoon pressure. Note however, that the overall 
normalisation of S v (I) depends sensitively on the value 
of the cocoon pressure. The effects of a higher pressure 
also include a steepening of S v (I) starting from the peak 
of this function towards the source core (Fig. [j]). This is 
not seen for variations of the viewing angle (Fig. In 
practice therefore, the best-fitting values for the cocoon 



Fig. 3. The influence of the viewing angle on the radio 
surface brightness. All model parameters as in Fig. [j] ex- 
cept for a variation in 9. Note that the measured length 
of the cocoon, Lj, changes from 86.6 kpc for the fiducial 
model to 98.5 kpc for 6 = 80° and 64.3 kpc for 9 = 40°. 




Fig. 4. Spectral index between 1.5 GHz and 5 GHz for the 
fiducial model and models with varying viewing angle, 9, 
and slope of the initial energy spectrum of the relativistic 
particles, p. All other model parameters as in Fig. [|. 

pressure is set mostly by the overall surface brightness of 
the cocoon while the viewing angle is mainly determined 
by the behaviour of S v (I) close to the core. 

Fig. H summarises the effects of the remaining model 
parameters on the surface brightness distribution. Chang- 
ing the aspect ratio of the cocoon, Rt, results in a scaling 
of S v (0 very similar to the effects of a changing viewing 
angle, 9. Decreasing ct2 leads to the end of the cocoon to 
become more blunt which in turn implies a lower value of 
Q!3, i.e. a slower backflow within the cocoon (Eqs. [L6| and 
|l6| ). This is analogous to an older particle population at 
a given distance from the hot spot and the resulting sur- 
face brightness distribution is similar to that of an older 
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Fig. 5. Effects of other model parameters on the radio sur- 
face brightness. Model parameters are as in Fig. |l| unless 
otherwise indicated in the legend. 

source. Changing the second shape parameter a\ has sim- 
ilar effects as that of a variation of oti and is therefore not 
shown. 

The main effects of changing the slope, p, and the high 
energy cut-off of the initial energy distribution of the rela- 
tivistic electrons, 7 max , is a change in the scaling of S„ (I) 
(see Fig. ||). This is caused by the dependence of the nor- 
malisation of the energy spectrum on p and 7max, Eq. 

Both scaling effects are similar to the effects of a 
variation of the viewing angle, 9. The effect of a steeper 
energy spectrum also causes an off-set in the distribution 
of the spectral index as a function of I as shown in Fig. |j. 
However, in the cocoon region closest to the source centre 
and least influenced by the hot spot emission, this effect 
is small compared to that of a variation of the viewing 
angle, 9. 

Finally, a different value of the slope of the gas density 
distribution in the source environment, /3, has a negligible 
effect on S v (I). The relevant curves for (3 — 1 and (3 = 2 
are indistinguishable from that for (3—1.5 shown in Figs. 

4-4- Degeneracy of parameters 

Both, ai and a2, control how pointed the shape of the co- 
coon is. The ends of the cocoon are dominated by the emis- 
sion from the hot spots which are not part of the model. 
Furthermore, if the jet direction is not stable over the life- 
time of the source, they may at times advance ahead of the 
rest of the cocoon which distorts the cocoon shape. This 
has been referred to as the 'dentist drill' effect (Scheuer 
1982). In the previous Sect, (see also Eq. [lfj) it was already 
shown that changing a\ and/or a.2 results in a change of 
the profile and magnitude of the backflow in the cocoon. 
This in turn causes changes to the surface brightness pro- 



files similar to changes of the source age. For these reasons 
the model cannot provide strong constraints on either of 
these parameters, at least not independent of the source 
age. Therefore I set in the following a,\ = 2 and a-i = 1/3. 

The surface brightness distribution predicted by the 
model is almost independent of the slope of the density 
profile of the external medium, (3. I therefore set (3 = 1.5 
without influencing the model results significantly. Note 
that because of Eq. ( |i"6| ) this then also implies a fixed 
value for 013 = 4.8. This may seem a rather significant 
restriction of the model but as was shown in the previous 
Sect, ai or ct2, through their influence on a^, and the age 
of the source, t, are degenerate model parameters. In the 
absence of geometrical constraints on ai or ct2, setting 
them to reasonable values allows the determination of the 
source age from the fitting method. The ages derived from 
the model will therefore always depend somewhat on the 
choice for the source geometry. 

The slope of the initial energy distribution of the rel- 
ativistic particles, p, mainly influences the spectral in- 
dex distribution in the cocoon. Observations suggest that 
2 < p < 3 and the model is therefore restricted to values 
in this range. It follows that the model is rather insensitive 
to the exact value of the high energy cut-off of the distri- 
bution, 7 m ax- A value of 7 max much smaller than the 10 4 5 
used in the fiducial model will, in addition to a change 
of the overall scaling, cause the emission region to shorten 
along the jet axis at a given observing frequency. However, 
because the highest frequency of the synchrotron radia- 
tion of the cocoon, f maXl depends on 7^ ax a small value of 
the high energy cut-off also implies a significantly smaller 
i/ max . Values substantially below i/ ma x ~ 100 GHz, which 
is assumed, here are unlikely in view of observations. I 
therefore set 7max = 10 4 5 . 

As was pointed out above, the effects of varying the 
initial slope of the energy distribution of the relativistic 
particles, p, are small in general. In Sect. || the 1 and 
2-dimensional comparison methods are applied to obser- 
vations of three FRII-type radio sources. It is shown there 
that in almost all cases the best-fitting model parameters 
for the 2-dimensional method require p to be close to 2. 
To prevent the degeneracy between p and 9 to influence 
the model fits in the 1-dimensional method which involves 
fewer degrees of freedom, I set p = 2 in this case. 

Finally, the degeneracy between i?x, the width of the 
cocoon, and the viewing angle, 9, is resolved by determin- 
ing the projected cocoon width R! T from the observed map 
and use Rt = R' T / sin# in the model calculations. 

Using these additional assumptions, the number of 
model parameters which are fitted by comparison to the 
observations decreases to four: The pressure in the cocoon, 
p c , the source age, t, the viewing angle, 9, and the initial 
slope of the energy distribution of the relativistic particles, 
p. In the case of the 1-dimensional comparison method p 
is fixed to a value of 2. 
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5. Application to Cygnus A, 3C 219 and 3C 215 

To test the model predictions for the source environment 
against direct X-ray observations over a range of differ- 
ent viewing angles, the radio data of three different FRII- 
type objects are used: the narrow-line radio galaxy Cygnus 
A, the broad-line radio galaxy 3C 219 and the radio-loud 
quasar 3C 215. According to orientation-based unification 
schemes of the various sub-classes of radio-loud AGN (e.g. 
Barthel 1989), the viewing angle, 9, of Cygnus A should 
be greater than those of 3C 219 and 3C 215. Furthermore, 
3C 219 and 3C 215 were selected because of their rather 
irregular radio lobe structure. The model is based on a 
very regular geometrical shape of the cocoon, Eq. (|l2|), 
and using 3C 219 and 3C 215 it is possible to estimate to 
what extent this restriction limits the applicability of the 
model. 

Note, that results of the model fitting for Cygnus A 
are presented in || and ^| only at one frequency. However, 
the model fits are always obtained for all sources using 
two maps at two different frequencies. 

5.1. Cygnus A 

For this source at z — 0.056 I used raw data from the 
VLA archive at 1.8 GHz and 5 GHz. The lower frequency 
observations were made in August/September 1987 in A- 
array and the 5 GHz data were obtained in January 1984 
in B-array. A detailed analysis of these and other obser- 
vations of Cygnus A can be found in Carilli et al. (1991). 
Standard calibration and self-calibration was performed 
using the software package AIPS. This resulted in two ra- 
dio maps with comparable angular resolution of 1.3" at 
the two observing frequencies. The CLEAN components 
were restored in maps with an individual pixel size of 
0.3" x 0.3". The map at 1.8 GHz with the hot spot emis- 
sion removed (see below) is shown in the upper panel of 
Fig. ^. The rms noise in the maps is 0.03 Jy beam -1 at 
1.8 GHz and 0.006 Jy beam -1 at 5 GHz. All pixels in the 
maps with a surface brightness below ha were discarded 
for the 2-dimensional comparison method. In the case of 
the 1-dimensional method only pixels below 3-er were ne- 
glected. The surface brightness distribution along the jet 
axis, S v (I), for the 1-dimensional method was obtained for 
both radio lobes along a cut from the core of the source 
to the eastern and western hot spot respectively. 

Since the hot spot emission is not modeled, it has to be 
removed from the maps and the 1-dimensional cuts. In the 
maps an aperture centered on the surface brightness peak 
in each lobe and with a radius of 2.6" corresponding to 
twice the beam width is removed (see Fig. ||). In the west- 
ern lobe the bright secondary hot spot (Carilli et al. 1991) 
and the bright ridge connecting this hot spot with the 
main one are also removed. For the 1-dimensional method, 
the distance of the hot spot to the edge of the lobe, AZ, 
was estimated as the distance of the maximum of S v (I) to 




Fig. 6. The 2-dimensional comparison method applied to 
Cygnus A. The upper panels shows the observed map at 
1.8 GHz. Contours are increasing linearly in steps of 0.43 
Jy beam -1 , starting at 0.14 Jy beam -1 , the 5-a noise level, 
to a peak just above 3 Jy beam -1 . The map is rotated by 
about 20° clockwise compared to the true position angle. 
The lower panel shows the x 2 -deviation of the best-fitting 
models. The filled contours show regions where x 2 = 
(white), < x 2 < 1, K X 2 < 2, 2 < x 2 < 10, 10 < x 2 < 
50, 50 < x 2 < 100 and x 2 > 100 (black). Note that the 
best-fitting model for the eastern and western lobes are 
not identical (see Table |J) . The results for both lobes are 
presented together for ease of comparison. The best-fitting 
model is obtained by fitting the 1.8 GHz map presented 
here and simultaneously the 5 GHz data (see Fig. 0) 



the last point at which this function has a value above 3a 
in the direction away from the source core. To remove the 
contribution of the hot spots to S v (I), all pixels within 
2AZ of the edge of the lobe were neglected in the following 
1-dimensional comparison process. 

The spatial resolution of the radio maps of Cygnus A 
is comparatively high. The 1.3" angular resolution corre- 
sponds to a spatial resolution of about 1.9 kpc. For many 
sources, particularly at high redshift, maps of such high 
quality are not available. In order to estimate the effects 
of a lower angular and therefore lower spatial resolution, 
I also convolved the two maps of Cygnus A with a Gaus- 
sian beam of 5" FWHM. In the case of the 2-dimensional 
method the radius of the aperture used to remove the 
hot spot emission was fixed to 5". For the 1-dimensional 
method, the surface brightness distribution along the jet 
axis was extracted from these lower resolution maps in 
the same way as in the higher resolution case. Note that 
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Fig. 7. Same as Fig. || but at 5 GHz. The contour levels 
in the upper panel are spaced linearly in steps of 0.12 Jy 
beam -1 starting at 0.03 Jy, the 5 — a level, to 0.76 Jy 
beam -1 . Higher contours are omitted for clarity. The \ 2 - 
contours are as in Fig. ||. 



the emission of the hot spots is smeared out over a larger 
area in the low resolution maps. To avoid any bias from 
the enhanced emission of the hot spot region I used the 
higher value of 2A/ obtained from the low resolution maps 
in both, the high and low resolution, 1-dimensional com- 
parison. 

5.1.1. The eastern lobe 

The eastern lobe of Cygnus A is covered by 46.6 indepen- 
dent telescope beams along the jet while in the widest part 
there are 22.2 beams across. In the lower resolution maps 
these numbers decrease to 12.9 and 6.1, respectively. Note, 
that only a fraction of the lobes has an observed surface 
brightness above the rms limits, i.e. they do not extend all 
the way from the hot spots to the core in the observations. 
This implies that the model fits are based on regions cov- 
ering less area than the theoretical extend of the cocoons. 
For both resolutions I find an axial ratio, R^, of 2.1 at 
the point where the lobe is widest from the 2-dimensional 
maps. The length of the lobe, Lj, is 64.4" for the low res- 
olution map and 60.6" for the high resolution map. The 
prediction of the best-fitting model in comparison with the 
observations is shown in the lower panels of Figs, bland 
for the 2-dimensional method and in Figs. || and g for 
the 1-dimensional method at 1.8 GHz. The 5 GHz data is 
not shown for the 1-dimensional method but is similar to 
the result at 1.8 GHz. The parameters of the best-fitting 
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Fig. 8. Comparison of the predicted surface brightness 
distribution along the cocoon with the low resolution ob- 
servations of the eastern lobe of Cygnus A. The crosses 
show the value of S v (V) at 1.8 GHz along the jet axis 
taken from the VLA map convolved with a 5" Gaussian 
beam. The solid line shows the best-fitting model with the 
model parameters given in table |[ Only data points left of 
the dashed line located at 2AZ from the tip of the cocoon 
(see text) were used in the fitting procedure to avoid the 
contribution of the hot spot. The entire observed range of 
S v (I) including the peak of the hot spot is shown in the 
inset. 

models are given in table 0. The errors on these and for all 
the following model fits are estimated using the boot-strap 
method (e.g. Press et al. 1992). It is not possible to cal- 
culate error estimates using the % 2 -values derived in the 
minimisation procedure directly as the values of the sur- 
face brightness in neighboring pixels are not independent. 
Roughly 2000 data sets were created by drawing data with 
replacement from the original set. The same minimisation 
procedure as in the original model fitting was then applied 
to them and the error given in the table is the la limit on 
the respective model parameters. 

Even after subtracting the contribution of the hot spot 
the deviations of the model from the observations are 
large. The fact that the model fit is poorer at 5 GHz is 
mainly caused by the smaller rms noise of the observed 5 
GHz map. As expected, the model cannot fit structures 
which appear as discrete surface brightness enhancements 
in the maps. This is particularly clear in the case of the 
bright arc seen just behind the hot spot in radio maps of 
the eastern lobe (see Figs. ^ and [?]) which also causes the 
secondary peak in S v (I) at about I — 0.83 (see Fig. ^|). 
Although convolving the maps with a larger beam 'draws' 
some flux from the hot spot into the arc mentioned above, 
the results for the two different resolutions are very simi- 
lar for both comparison methods. At both frequencies the 
observed maps show a concentration of the radio emission 
towards the centres of both lobes. This region is also ex- 
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Table 3. The best-fitting model parameters from the 1 and 2-dimensional comparison methods. 
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Fig. 9. Comparison of the predicted surface brightness 
distribution along the cocoon with the high resolution ob- 
servations of the eastern lobe of Cygnus A. The data is 
taken from the unconvolved VLA map. As in Fig. [§] the 
solid line shows the best-fitting model and only data left 
of the dashed line was used in the fitting procedure. 



tended a long way along the jet axis, particularly in the 
western lobe (see following section). In Figs. || and ^ it is 
clear that the model cannot fit this concentration prop- 
erly. The very smooth lobes of the model are 'fatter' than 
the observed lobes further away from the hot spots and do 
not extend as far back as the observations indicate. This 
is particularly striking in the western lobe at 5 GHz (see 
Fig. (?]). This clearly illustrates the limitations of the model 



assumption of a smooth shape of the cocoon and a regular 
backflow within the cocoon. The large uncertainty of the 
viewing angle, 9, is caused by the model mainly depend- 
ing on sin 9 which changes only by a factor 1.2 within the 
estimated errors. This is comparable to the uncertainties 
of the other model parameters. 

Using Eqs. (||) and (||) The power of the jet, Q a , and 
the central value of the density distribution of the gas 
surrounding Cygnus A, p , are calculated. The results are 
given in Table |j. Here I assume that the core radius of 
the environmental density distribution in Eq. (||) is given 
by a = 24 kpc. The viewing angle 9 is inferred from the 
flux ratio of the jet to the counter-jet in the two lobes of 
Cygnus A (Hardcastle et al. 1999). This assumes that the 
two jets are identical and that the observed flux ratio is 
entirely due to relativistic beaming effects. Furthermore, 
a constant bulk velocity within the jets, i>j, is assumed 
and set to 0.62 c. Variations of Vj across the source and 
asymmetries between the two jet sides will significantly 
influence the estimate for 9. However, the model is consis- 
tent with the estimate given by Hardcastle et al. (1999). 

The observed central density of source environment 
given in Table |] is derived from X-ray observations 
ROSAT of the hot gas surrounding Cygnus A Hardcas- 
tle & Worrall (2000). For this, the prescription of Birkin- 
shaw & Worrall (1993) for the conversion of central sur- 
face brightness to central proton density was used. The 
core radius, a Q was estimated by Hardcastle & Worrall 
(2000) from the X-ray observations and I use their value, 
a Q = 24 kpc, in converting from the density parameter 
Po a S given by the model to p . The value thus found from 
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Table 4. Properties of the source environment derived from the best-fitting model parameters in comparison with 
observations. For the determination of p the core radii of Hardcastle & Worrall (2000) inferred from X-ray observations 
were used. These are a Q = 24 kpc for Cygnus A, a a — 140 kpc for 3C 219 and a Q = 204 kpc for 3C 215. 
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the best-fitting model agrees within the error with the X- 
ray observations. 

5.1.2. The western lobe 

The western lobe of Cygnus A is covered by 53.5 indepen- 
dent beams along the jet and 23.2 beams at the widest 
point perpendicular to the jet. For the lower resolution 
maps I find 14.6 and 6.3 beams, respectively. Similar to 
the eastern lobe the observed emission does not extend all 
the way from the hot spot to the core and so the model fit 
is based on a smaller area. From the 2-dimensional maps I 
find R' T = 2.3 for the widest part of the lobe and the length 
of the lobe is 73.0" and 69.5" for the low and high resolu- 
tion case respectively. The best-fitting models for the two 
different resolutions again agree well. The models yield an 
age for the western lobe somewhat higher than that of its 
eastern counterpart (see Table which is mainly caused 
by its greater length. However, the pressure within the co- 
coon is remarkably similar in both lobes. This implies also 
good agreement between the estimates for the jet power 
and the density of the source environment between the 
two sides of Cygnus A (Table f|) . 



5.2. 3C 219 

For 3C 219 at z = 0.1744 VLA maps at 1.5 GHz in B-array 
and 4.9 GHz in C-array were used. The observations were 
taken in October and December 1998 by Dennett-Thorpe 
et al. (in preparation) who also performed standard re- 
duction on the data set using AIPS. The resolution of 
the resulting maps is roughly 4.3" and the rms noise is 
2.9 ■ 10~ 4 Jy beam" 1 at 1.5 GHz and 4.4 ■ 10~ 5 Jansky 
beam -1 at 4.9 GHz. The northern lobe of 3C 219 has a 
rather irregular shape and no clear hot spot (Clarke et al. 
1992). This leads to large ambiguities in the determination 
of its length or the geometrical parameters needed for the 
model presented here. I therefore only used the southern 
lobe which has a length of 40.4" and an aspect ratio at 
its widest point of R' T — 1.6. The lobe is covered by 9.4 
independent beams along the jet axis and by 5.9 beams 
perpendicular to it. The jet and counter-jet in 3C 219 are 
unusually bright and so the jet emission was removed from 
the maps of the southern lobe. For the 1-dimensional com- 
parison method, I extracted S u (I) along a line off-set by 
2" to the south of the line connecting the core of the source 
with the hot spot of the southern lobe. This avoids con- 



C. R. Kaiser: Environments and ages of radio sources 



15 



tamination of the surface brightness distribution by the jet 
emission. The hot spot in the southern lobe of 3C 219 is 
somewhat set back from the edge of the lobe. An aperture 
with a radius of 4.3", i.e. the width of the telescope beam, 
centered on the surface brightness peak was removed from 
the maps. For the 1-dimensional comparison, only values 
of S v (I) core-wards of the hot spot were used. 

The parameters of the best-fitting model are given in 
Table |[ The uncertainties of the model parameters is com- 
parable to those found for the two lobes of Cygnus A. The 
angle to the LOS of 3C 219 is found to be smaller than 
that of Cygnus A. Since 3C 219 is a broad line radio galaxy 
this is in the expected sense, but the value of 9 predicted 
by the model is about double that inferred from the flux 
ratio of the jet and counter-jet. As mentioned above, the 
jets of 3C 219 are unusually bright and this may reflect 
some enhanced disruption of the jet flow by turbulence or 
even a complete restart of the jets in this source (Clarke & 
Burns 1991). The latter possibility has let Schoenmakers 
et al. (2000) to include this source among their examples of 
Double-Double Radio Galaxies (DDRG). The morphology 
of these sources strongly suggests restarting jets (Kaiser et 
al. 2000). The large jet to counter-jet flux ratio of 3C 219 
may therefore be caused by effects other than relativistic 
beaming. The best-fitting value of 9 = 65° is consistent 
with orientation-based unification schemes. 

The central density of the gas surrounding 3C 219, p Q , 
predicted by the model is considerably smaller than that 
inferred from X-ray observations. To derive p I used a a 
I I kpc (Hardcastle & Worrall 1999). This discrepancy 
will be discussed in Sect. ^| 

5.3. 3C 215 

This source is a radio- loud quasar at z = 0.411 with very 
irregular morphology (Bridle et al. 1994). I obtained raw 
observational data from the VLA archive at 1.5 GHz in 
A-array and 4.9 GHz in B-array. The 1.5 GHz observa- 
tions were carried out by Miley in May 1986 while the 
4.9 GHz observations were taken by Hough in December 
1987. Again standard reduction with AIPS was performed 
on the data and resulted in two maps with an angular 
resolution of 1.9". The maps were restored using a pixel 
size of 0.3" xO. 3" and the rms noise is 1.6 ■ 10 -4 Jansky 
beam -1 at 1.5 GHz and 4.0 • 10 -5 Jansky beam -1 at 4.9 
GHz. The southern half of 3C 215 is very distorted with 
the jet bending in various places with enhanced surface 
brightness (Bridle et al. 1994). This part of the source is 
not consistent with a regular FRII-type lobe morphology 
and resembles in its outer regions an FRI-type structure. 
Therefore no attempt was made to apply the model to the 
southern part of the source. The northern lobe is more 
regular, however, the hot spot here is weak and the lobe 
widens considerably close to the core in a north-eastern 
direction. The lobe has a length of 26.6" and its aspect 
ratio R' T at the point where the width of the lobe is great- 



est is 1.2. The hot spot in the northern lobe is not located 
at the very edge of the lobe similar to the southern lobe 
of 3C 219. An aperture centered on the surface bright- 
ness peak with a radius 3.8" corresponding to the size of 
two telescope beams was removed from the map. For the 
1-dimensional comparison the surface brightness distribu- 
tion was extracted along a line off-set by about 2" to the 
east from the core-hot spot direction to avoid emission 
from the jet. Again only pixels core-wards from the hot 
spot were used in the 1-dimensional case. The northern 
lobe of 3C 215 is covered by 14.0 independent beams along 
the jet and 11.7 beams at the widest point perpendicular 
to the jet. 

Parameters of the best-fitting model and the derived 
properties of the environment of 3C 215 are given in Ta- 
bles H and ||. The uncertainties for the model parameters 
are considerably larger for this source than for the two pre- 
vious ones. The viewing angle to the jet axis, 9, is smaller 
than for Cygnus A or or 3C 219. This is again consistent 
with the predictions of unification schemes as 3C 215 is a 
quasar. The smaller observed value is again inferred from 
the flux ratio of the jet and counter-jet of 3C 215 (Bri- 
dle et al. 1994). Similarly to 3C 219 this ratio may be 
increased in 3C 215 because of the distorted morphology 
of its large scale radio structure. The southern jet does 
not seem to be embedded in a cocoon and it is therefore 
unlikely that the two jets are intrinsically identical. The 
flux ratio probably reflects physical processes other than 
purely relativistic beaming. 

For the determination of p I used a = 204 kpc from 
Hardcastle & Worrall (1999). The density of the gaseous 
environment of 3C 215 is predicted to be much lower than 
inferred from X-ray observations. Discussion of this point 
is deferred to Sect. ^. 

6. Discussion 

6.1. 1-D versus 2-D method 

In general it is expected that an increase in the amount 
of information available to a given model fitting proce- 
dure should decrease the uncertainty in the determina- 
tion of the model parameters. Comparing the results of 
the 1 and 2-dimensional methods for the lobes of Cygnus 
A one would therefore expect that the parameter uncer- 
tainties decrease for the higher resolution maps. This is 
indeed the case for the 2-dimensional method, particularly 
for the western lobe. However, the 1-dimensional method 
shows the opposite. Here the uncertainties are larger for 
the higher resolution maps. This is caused by the attempt 
to fit a very smooth model for the surface brightness to 
observational data which shows considerably greater lo- 
cal variations than the model. The fact that the model 
is unable to fit the local surface brightness structure ob- 
served was already noted in the previous Sect.. In the 2- 
dimensional case the off-axis parts of the cocoons provide 
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additional information and the influence of local structure 
is therefore to some extent averaged out in the fitting pro- 
cedure. In other words, the 2-dimensional method is able 
to make use of the larger amount of information in higher 
resolution maps. The 1-dimensional method is restricted 
to a cut through the lobe. Here, the averaging effect of 
a larger telescope beam provides for a smoother surface 
density profile and the model fits the data better. 

In the case of radio maps of low resolution which do not 
or only barely resolve the lobes perpendicular to the jets 
along most of their lengths, the 2-dimensional maps add 
little information to 1-dimensional cuts along the lobes. 
In these cases the additional model parameter p in the 
2-dimensional method presented here allows the model to 
fit the data with a large range of parameter combinations. 
The uncertainties of the model parameters are then larger 
than for the 1-dimensional method. This effect can be seen 
for 3C 215 and, to a lesser extent, for the low resolution 
maps of the western lobe of Cygnus A. It is of course pos- 
sible to fix the value of p for the 2-dimensional method as 
well but this does not significantly improve the constraints 
on the model parameters compared to the 1-dimensional 
method. 

Which method is the best to use for a particular set of 
radio maps depends on the quality of the maps. If the ra- 
dio lobes are well resolved along to the jet axis as well 
as perpendicular to it then the 2-dimensional methods 
will provide better constraints. However, it is computa- 
tionally expensive. In the case of poorer resolution the 2- 
dimensional method will not add anything to the results 
obtained from a 1-dimensional comparison. For heavily 
distorted lobe structures both methods will fail but the 1- 
dimensional method may still provide order of magnitude 
estimates if the lobes are not entirely dominated by bright 
localised structure. 

6.2. Determination of viewing angles 

Orientation-based unification schemes attempt to explain 
radio galaxies and radio-loud quasars as essentially the 
same type of objects albeit viewed at different angles to 
the jet axis (Barthel 1989). The broad line radio galaxies 
may then be identified as the low redshift analoga to radio- 
loud quasars. These unification schemes imply that the 
viewing angle, 8, to the jet axis is greater than about 45° 
for radio galaxies and smaller for quasars. Tests of this 
hypothesis are often inconclusive as the determination of 
8 is difficult in practice. 

As described above, the viewing angle 8 may be deter- 
mined from radio observations alone. However, since the 
model depends mostly on sin6>, the uncertainties are large 
particularly for 8 > 45° . Fortunately, the other model pa- 
rameters and the source and environment properties in- 
ferred from these do not depend strongly on 8 and its 
error in this range. This agrees with the results of Wan & 
Daly (1998) who study the effects of the viewing angle on 



the determination of a variety of source properties in great 
detail. The source orientation can only be determined ac- 
curately if the 2-dimensional comparison method is used 
on radio maps well resolved perpendicular to the jet axis. 
This agrees with the findings of Wan & Daly (1998) who 
use a different source model. Despite the large uncertain- 
ties in the model results presented here it is interesting 
to note that for the three sources studied the best-fitting 
values of 8 are smaller for the quasar (3C 215) and the 
broad line radio galaxy (3C 219) than for the radio galaxy 
(Cygnus A). 

6.3. Environments of FRII sources 

As I showed above, the model parameters of the best- 
fitting models can be used to infer the density parame- 
ter, p a^, describing the density distribution in the envi- 
ronment of FRII objects. It is impossible using just this 
model to separate this into the central density, p a , and 
the core radius, a Q . Furthermore, the model is insensitive 
to the slope, f3, of the external density profile (see Sect. 
4.3). To infer central densities, values for the core radius 
and (3 must be taken from other observations. In Table 
^ p for the three example sources is shown for the core 
radii given by Hardcastle & Worrall (1999) and (3 — 1.5. 
As noted in Sect. 0, the central densities found for 3C 219 
and 3C 215 are inconsistent with those derived from the 
X-ray observations. 

In several studies it was found that the pressure of the 
FRII source environment derived from X-ray observations 
apparently exceeds the pressure inside the radio lobes (e.g. 
Hardcastle & Worrall 2000 and references therein). The 
discrepancy in density found here is essentially the same 
phenomenon in the framework of isothermal density dis- 
tributions for the source environments as described by the 
(3- model, Eq. (Q). All recent models for the evolution FRII 
sources, including the one discussed here, are based on the 
assumption that the cocoon is overpressured with respect 
to its surroundings. The X-ray observations seem to con- 
tradict this assumption. 

Can the discrepancies be resolved? In the model de- 
scribed here I approximate the density distribution in the 
source environment, assumed to follow the /3-model, by 
power laws. For cocoons extending well beyond the core 
radius the exponent of the power law is given by (3 = 3(3'. 
In the case of Cygnus A this is justified as both lobes ex- 
tend over more than 100 kpc while a — 24 kpc. In the 
analysis f3 = 1.5 was used which is identical with the result 
of Hardcastle, (3' = 0.5. For 3C 219 and 3C 215 (3' = 0.9 
from the X-ray observations and in both sources the ra- 
dio lobes only extend to about 1.3a . Formally the model 
cannot be applied to these two sources because the under- 
lying dynamical model is valid only for (3 < 2 or /3' < 2/3. 
However, as the cocoons in both sources do not extend 
much beyond the core radius, the density distribution in 
the environments of 3C 219 and 3C 215 are well repre- 
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sented by a power law with exponent (3 = 1.5. In these 
cases po found from the model is not identical to the cen- 
tral density of the /3-model. From Eq. p) it can be seen 

3/3' /2 

that the latter is given by p (l + x 2 ) = 3.8p D , where 
x is the length of the radio lobe in units of a a . This cor- 
rection is insufficient to make the model results consistent 
with the X-ray observations of 3C 219 and 3C 215. 

X-ray observations of the hot gaseous environment of 
AGN are influenced by the presence of the active galaxy. 
Before the environment properties can be extracted, the 
bright X-ray emission of the AGN itself appearing as a 
point source must be carefully removed. In the case of 
radio-loud objects the large scale structure caused by the 
jets can also alter the X-ray emission of its surround- 
ings. The magnitude of these effects is difficult to esti- 
mate if the X-ray observations do not fully resolve the 
scale of the cocoon. The hot spots at the end of the jets 
are strong sources of inverse Compton scattered X-ray 
photons. When resolved, these inverse Compton scattered 
photons are found to distort the X-ray contours of the 
extended emission (Cygnus A, Carilli et al. 1994; 3C 295, 
Harris et al. 2000). The more extended cocoon material it- 
self may also act as a scatterer of CMB photons or of AGN 
emission (Brunetti et al. 1999). Finally, the bow shock 
surrounding the cocoon compresses and heats the gas in 
the source environment. Kaiser & Alexander (1999b) give 
an estimate for the expected X-ray luminosity from this 
shocked layer of gas, 
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where two typographical errors are corrected. Expressions 
for I 2 /a2, ai and P/R may be found in Kaiser & Alexan- 
der (1999b). The square brackets mean the difference of 
the exponential function at the limits of the observing 
band v\ and v% in the rest frame of the source. The re- 
sulting X-ray luminosity can be converted to a ROSAT 
count rate using the internet version of PIMMSQ. For 
the best-fitting model parameters for Cygnus A this im- 
plies that roughly 2% of all counts of the ROSAT ob- 
servations presented by Hardcastle & Worrall (2000) at- 
tributed to the entire extended emission come from the 
layer of shocked gas in between bow shock and cocoon of 
this source. For 3C 219 this region around the southern 
lobe alone contributes about 10% of the relevant counts. 
The contribution of the northern lobe of 3C 215 is negli- 
gible (~ 0.1%). However, in this case at least three very 
compact radio emission regions along the distorted jets are 
detected which may be powerful inverse Compton sources 

1 PIMMS was programmed by K. Mukai at the High En- 
ergy Astrop hysics Science Archive Research Center of NASA, 
available at http://heasarc.gsfc. nasa.gov /Tools /w3pimms.htm\ 



(Bridle et al. 1994). Only X-ray observations resolving the 
scale of the radio structure in this and other sources will 
allow us to determine the contribution to the total X-ray 
emission from such compact regions and hot spots. 

Even without detailed observations it is clear that pow- 
erful radio sources can contribute significantly to the ex- 
tended X-ray emission in the central part of their envi- 
ronments. Helsdon & Ponman (2000) show for the case 
of loose groups of galaxies that such an overestimate at 
the centre of an X-ray surface brightness profile may lead 
to overestimates of the core radius and also of (3' when 
fitted with a /3-model. Using these overestimated values 
of a a and (3 then yields values for p which are too low 
compared to the 'real' central density. An exact analysis 
of the magnitude of this effect is beyond the scope of this 
paper. In any case, this effect may explain the discrep- 
ancies found between the two methods of estimating the 
central densities and pressures of the FRII environments. 

7. Conclusions 

In this paper an analytical model for the surface bright- 
ness distribution of the cocoons of FRII-type radio sources 
is developed. The model is based on the self-similar model 
of Kaiser & Alexander (1997, KA) and its extension by 
Kaiser et al. (1997, KDA). The cocoon is split into small 
volume elements, the temporal evolution of which is traced 
individually. These elements are then identified with cylin- 
drical slices which are rotationally symmetric about the 
jet axis. The bulk backflow of these slices is determined 
self-consistently from the constraints on the cocoon shape. 
Thus a 3-dimensional model of the synchrotron emissivity 
is constructed. Projecting this along the LOS, the model 
yields surface brightness profiles in one or two dimensions 
in dependence of several model parameters. It is shown 
that a number of degeneracies prevent constraints to be 
placed on all parameters. Comparatively robust estimates 
can be found for the cocoon pressure, p c , the source age, t, 
the angle of the jet to the LOS, 9, and, if the 2-dimensional 
comparison method is used, the initial slope of the energy 
distribution of the relativistic particles, p. 

The model may be viewed as an extension of the clas- 
sical spectral aging methods. However, it accurately takes 
into account the loss history of the relativistic particles 
and the evolution of the magnetic field in the cocoon. The 
age estimates derived from the model are therefore accu- 
rate and not only lower limits. It is also shown that diffu- 
sion of relativistic particles cannot significantly distort the 
energy spectrum of the relativistic particles in the cocoon. 

In Sect. H the model predictions were compared with 
observational data for the three sources Cygnus A, 3C 
219 and 3C 215. The four free model parameters are con- 
strained to within at least a factor 2. The uncertainties of 
the cocoon pressure are remarkably small, less than 50%, 
while the source ages are constrained by the model to 
within 80% or better. The errors on the viewing angle are 
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large as the model mainly depends on sin 9 which changes 
in the range of orientation angles found here (31° to 90°) 
by only 52%. The uncertainties in the model parameters 
translate to errors in the determination of the jet power, 
Qo, and the density parameter, p »f • However, the result- 
ing uncertainties still allow at least order of magnitude 
estimates for these quantities. 

The model results do not depend strongly on the res- 
olution of the radio observations used in the fitting pro- 
cedure. The results for both lobes of Cygnus A derived 
at 1.3" and 5" resolution are consistent with each other 
within the errors. However, the uncertainties in the param- 
eter determination increase significantly with decreasing 
resolution. For radio maps which barely resolve the lobes 
perpendicular to the jet axis, the 2-dimensional method 
adds little in terms of constraints on model parameters 
to the 1-dimensional method. From the examples stud- 
ied here it is not possible to give a firm lower limit for 
the resolution of maps used to fit with the model. How- 
ever, the parameter uncertainties calculated here suggest 
that at least 4 independent beams alo ng the jet axis are 
required. As was pointed out in Sect. 3A, distorted lobe 
structures may indicate large turbulent flows within the 
cocoon. As the model is based on the assumption of largely 
non-turbulent backflow it cannot be used for such irregu- 
lar sources. 

The model can be used to infer the viewing angle of 



FRII-type sources to the LOS. Despite the large uncer- 



some of the questions related to the cosmological evolu- 
tion of the FRII population and unification schemes of 
radio-loud AGN. 
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Appendix A: Diffusion of relativistic particles in 
the cocoon 

A.l. Quasi-linear anomalous diffusion 

In the model it will be assumed that the magnetic field 
within the cocoon is completely tangled on a scale r t much 



smaller than the size of the cocoon, L- y In this case the 
relativistic electrons can diffuse quickly within 'patches' 
of coherent magnetic field of a size r t . However, it is diffi- 
cult for them to escape their patch as this would involve 
diffusion perpendicular to field lines. In the presence of 
turbulent motion within the plasma underlying the mag- 
netic field, patches can be stretched out and in this case 
the diffusion of charged particles into an adjacent patch 
becomes more likely. Duffy et al. (1995) calculate how far 
a given electrons has to travel along a field line in its orig- 
inal patch before the patch becomes so stretched out that 
it can escape. Note here that this treatment is valid only 
in the quasi-linear regime for which the relative ampli- 
tude of irregularities of the magnetic field is much smaller 
than the ratio of the turbulent correlation lengths perpen- 
dicular and parallel to the local magnetic field. The effi- 
ciency of this diffusion then depends crucially on whether 
the electron is able to travel this distance ballistically or 
whether it must diffuse along the field line. The more ef- 
ficient ballistic regime requires that 

r g > V«L*e) (A.l) 

where r g is the gyro-radius of the particle in the field, K± 
is the diffusion coefficient perpendicular to the field and 
t e is the time it takes the particle to escape the patch. For 
a relativistic electron moving at speed v and correspond- 
ing Lorentz factor 7 following Chuvilgin & Ptuskin (1993) 
K± = £Kb/ (e + 1) with kb = -fv 2 m c c/ (3eB) the Bohm 
diffusion coefficient and e = v CQ \jv % , where v co \ is the rate 
of collision of the particle and v g is its gyro-frequency. For 
the assumption of Blundell & Rawlings (2000) that the 
coherence length of the magneti c fie ld is roughly 10 kpc 
in all directions, the inequality ( A.l ) yields e < 7 • 10~ 12 
for a magnetic field strength of 130 fiG, appropriate for 
the fiducial model of Sect. 4.5, and an electron with a 



Lorentz factor 1000. This implies that the time between 
collisions of this electron must be greater than 2 Myr, i.e. 
the mean free path of the electron is of order 600 kpc. 
This is clearly unphysical in the case of Lj ~ 100 kpc as 
relativistic particles would then simply escape the cocoon 
very quickly. 

The relativistic electrons in the cocoon must therefore 
diffuse along the magnetic field lines in between jumps 
from one patch of coherent magnetic field to another. The 
expression of Duffy et al. (1995) for the rms diffusion 
length after a time t, x, in this case can be approximated 
by 



r„ct 



3elog(V2e) 



(A.2) 



Assuming that the mean free path of the electron is less 
than the coherence length of the magnetic field, i.e. less 
than 10 kpc, I find for the same magnetic field and Lorentz 
factor x<3kpcifi = l Myr. A significant mixing of 
relativistic particles along the cocoon due to anomalous 
diffusion is therefore unlikely. 
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A. 2. Non-linear diffusion 

Of course, it may be argued that as the above analysis only 
applies to the quasi-linear regime, the diffusion timescale 
in a highly turbulent flow may be much shorter. Consider 
such a flow to be present in the cocoons of FRII-type 
radio sources. In this case, the relativistic electrons may 
diffuse through the observed lobes within a time short 
compared to the age of the source. However, since diffusion 
is a stochastic process and the geometry of the cocoon 
is elongated, most of the particles will leave the cocoon 
sideways and will be lost to the surrounding gas before 
traveling large distances along the cocoon. It is likely that 
the diffusion time for the relativistic particles depends on 
their energy and so the diffusion losses, if present, will 
significantly change the emission spectrum which is not 
observed (e.g. Roland et al. 1990). Even in the case of 
efficient energy independent diffusion the observed radio 
lobes should show diffuse edges in low frequency radio 
maps. Again this is not observed (e.g. Roland et al. 1990, 
Blundell et al. 2000a,b). 

From the above it is clear that in the presence of very 
efficient diffusion some special confinement mechanism for 
the relativistic particles in the cocoon preventing their es- 
cape sideways is needed. This may be provided for by the 
compression and shearing of the tangled magnetic field at 
the edges of the cocoon. This process will align the mag- 
netic field close to the cocoon edge with this surface and 
therefore act as a kind of magnetic bottle. The order thus 
introduced in the originally tangled magnetic field due to 
this process leads to an enhanced polarisation of the emit- 
ted radiation in this region (e.g. Laing 1980). 

Large volume compression ratios are ruled out as the 
sound speed in the cocoon is high (e.g. KA). However, for 
a conservative estimate consider a volume compression ra- 
tio of 10 in the sheet of compressed material at the edge 
of the cocoon of a given source. This already implies that 
the maximum theoretical value of polarisation of the syn- 
chrotron emission should be observed at the edge of the 
radio lobes (Hughes & Miller 1991). For adiabatic com- 
pression of a tangled magnetic field the strength of the 
field increases by a factor of roughly 4.6. The rate of col- 
lisions of a given particle is probably increased as well 
as the irregularities in the magnetic field are also com- 
pressed. However, the case of v co \ = const, provides for a 
lower limit of e and so with k± = €Kb/ (e + 1) it is clear 
that the diffusion coefficient within the compressed region 
perpendicular to the field lines decreases by a factor 25 
at most compared to the inner cocoon. This implies that 
in a given time a relativistic particle diffuses a 5 times 
shorter distance in the compressed boundary layer com- 
pared to the inner regions of the cocoon. To significantly 
influence the time the particle remains in the cocoon this 
compressed layer must at the very least occupy 20% of 
the width of the cocoon which is unlikely. Note, that al- 
though the magnetic field is compressed in the boundary 



layer and therefore aligned with the cocoon surface some 
field lines may still be perpendicular to this surface. Along 
these field lines diffusion will be even faster allowing many 
electrons to escape the cocoon. 

From this I conclude that diffusion, even if it is highly 
effective in the inner cocoon and there exists a compressed 
boundary layer along the edges of the cocoon, will not 
alter the distribution of relativistic particles within the 
cocoon. This allows us to use the spatial distribution of the 
synchrotron radio emission of FRII sources to infer their 
age. The model developed in the following can be viewed 
as an extension to the classical spectral aging methods in 
that it takes into account the evolution of the magnetic 
field in the lobe 



